Coherent atom-trimer conversion in a repulsive Bose-Einstein condensate 



O 

o 

(N 



(N 



> 

m 
O 
i> 
o 

^ : 
Oh: 

CT 1 - 



13 



H. Jing 1 , J. Cheng 2 , and P. Mcystrc* 1 

1 Department of Physics, The University of Arizona, Tucson, Arizona 85721, USA 
2 School of Physical Science and Technology, South China University 
of Technology, Guangzhou 510640, People's Republic of China 
(Dated: February 1, 2008) 

We show that the use of a generalized atom-molecule dark state permits the enhanced coher- 
ent creation of triatomic molecules in a repulsive atomic Bose-Einstein condensate, with further 
enhancement being possible in the case of heteronuclear trimers via the constructive interference 
between two chemical reaction channels. 
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The experimental realization of Bose-Einstein conden- 
sates (BEC) in dilute atomic vapors has led to a num- 
ber of spectacular advances with implications well past 
the confines of traditional atomic, molecular and opti- 
cal (AMO) physics [1], and with profound relevance for 
fields ranging from condensed matter physics to quantum 
information science. A development of particular rele- 
vance for the present study is the generation and probing 
of ultracold molecular samples and of molecular conden- 
sates from atomic BEC, using magnetic Feshbach reso- 
nances (FR) and photoassociation (PA) [2-6] to control 
the dynamics of the system. In particular, it has recently 
been shown that these two techniques can be combined 
to achieve the efficient and stable conversion of atoms 
into molecular dimers [3]. 

The goal of this paper is to demonstrate theoretically 
that these techniques can in principle be extended to the 
generation of molecular trimers. An important new re- 
sult is that the creation of heteronuclear trimers can be 
significantly enhanced by the constructive interference of 
two quantum channels leading to their formation from a 
two-component atomic condensate. The basic idea is to 
first create highly excited dimers via a standard FR, and 
then to couple them to a bound trimer via PA. A key 
element of the scheme is to exploit a coherent population 
trapping (CPT) technique to prevent the dimer popula- 
tion from becoming significant throughout the conversion 
process. Such a scheme has previously been proposed for 
the creation of molecular dimers [3] , and has been theo- 
retically demonstrated to be stable for a broad range of 
conditions. 

This proposed extension of "superchemistry" [6] ma- 
nipulations from dimers to trimers exploits the existence 
of three-body bound states in ultracold atomic samples 
as the scattering length for two-body collisions becomes 
infinite. As discussed in Ref. [7] this occurs not only 
for identical particles, a situation first considered by Efi- 
mov [8-10], but also for two identical and one different 
particle. In particular, in the case of two heavy particles 
and one light particle these authors found that the heavy- 
heavy interaction does not matter, and there is an infinite 
number of bound trimer states as long as the light-heavy 



scattering length becomes infinite. In the case of the 
interaction between one heavy and two light particles, 
by contrast, there is no infinite number of bound states 
unless both two-body scattering lengths become simulta- 
neously infinite. 

We mentioned that a point of particular interest in 
the formation of heteronuclear trimers A2B is the role 
of quantum interferences between the two paths that in- 
volve intermediate dimers A2 (path AA) and AB (path 
AB), respectively. To set the stage we consider first the 
path A A, whereby A A dimers formed via a FR are pho- 
toassociated with an atom B to form the trimer A2B [3]. 
The essence of the idea is to minimize the occupation 
of intermediate dimers by exploiting an atom-molecule 
dark state that permits the direct association of atoms 
into trimers without the formation of a substantial dimer 
population in the process [5]. We note that while the 
proposed scheme is experimentally challenging, recent 
progress in the manipulation of dimer-atom resonances 
[11-14] indicates that it might become realizable in the 
near future. Note also that this process is quite different 
from a FR-induced dimer-trimer mixture in a resonant 
condensate [15]. 

Our model system consists of a Bose-Einstein conden- 
sate of atoms coupled to molecular dimers via a FR, these 
dimers being in turn photoassociated to the atoms to 
form bound trimers. Denoting the strength of the atom- 
dimer coupling with detuning 5 by , the Rabi frequency 
of the PA laser by and its detuning by A, the dynam- 
ics of the system is described at the simplest level by the 
model Hamiltonian 
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Here, the annihilation operators tpi(r), where the indices 
i,j = a,b,d,g stand for atoms (A and B), dimers and 
trimers, satisfy bosonic commutation relations, and the 
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collision terms proportional to Xij describe s-wave colli- 
sions between these species. (Trimer formation via the 
path AB only can be investigated in a similar fashion.) 
The nonlinear bound-bound coupling between dimers 
and trimers is typically induced by a narrow-frequency, 
continuous-wave PA laser, for which the Franck-Condon 
factor can be calculated by resonant scattering theory [2, 
11] and may be tuned by pumping methods [16-17]. 

We assume in the following that the main features of 
the dynamics are adequately described by a mean-held 
analysis, ipi — > ^/nipi, where n is the initial atomic den- 
sity. In this limit, the system is described by the equa- 
tions of motion (h = 1) 
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where Ai = \[y/n, f2i = Q^y^n and the decay rate 7 ac- 
counts for the loss of untrapped dimers. To reduce these 
losses we exploit a CPT technique that relies on the exis- 
tence of an approximate atom-molecule dark state. Such 
techniques arc well known in the case of linear systems, 
where they permit the transfer of population from an ini- 
tial to a final state via an intermediate state that remains 
unpopulated at all times. This is the basis for stimulated 
Raman adiabatic passage (STIRAP), which achieves this 
goal via a so-called counter-intuitive sequence of pulses 
[2]- 

CPT and STIRAP rely explicitly on the validity of the 
adiabatic theorem, which applies only to linear systems. 
Unfortunately, in the situation at hand two-body colli- 
sions render the system nonlinear, and it is not immedi- 
ately obvious that STIRAP still works in that case. This 
problem has recently been investigated in Ref. [4], which 
shows that an approximate adiabatic condition can still 
be achieved by linearizing the nonlinear system around 
the intended adiabatic evolution. If the cigcnfrcqucncies 
of the linearized system are real, the associated "normal 
modes" will not grow in time, and the system is stable. 
Hence a system initially prepared in a CPT state will ap- 
proximately remain in that state at all times, although 
the adiabaticity condition may be difficult to fulfill at the 
later stages of the evolution. In contrast, if the eigenval- 
ues of the linearized problem are complex, the system is 
dynamically unstable for some parameter values [3] under 
which adiabaticity breaks down. 

For the specific case of trimer formation it is easily 



shown that under the generalized two-photon resonance 
condition 

A = -5 + 2(2Xag + Xbg ~ Xgg)nN g , s 

+ (4Xaa - 2Xag + ^Xab + Xbb ~ Xbg)nN a ,s (2) 

cqs. (fTJ admit a steady-state solution with no dimer pop- 
ulation, 
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where i = 1 and k = 4. (For the path AB, Eq. (3) 
remains the same, but with i = 2 and k = 1, so that 
A^B < ArAA for the 

same external parameters. Note 1 
that Nd, s = and the same two-photon resonance con- 
dition holds in both cases.) 

This suggests that approximate CPT dynamics such 
that the dimer population Nd remains small at all times 
can be achieved for an appropriate "counter-intuitive" 
time dependence of the laser detuning A [3] (see also 
Ref. [18]). Additionally, it may be possible to exploit a 
feedback technique [19] to further stabilize this process. 

0.8 



c 0.6 

Q 

a 0.4 

o 

CL 



0.2 



0.8 



e 0.6 

o 

TO 

3 _ . 



0.2 







path (AA) 


atom A (b) 


CPT 


atom A (a) ~~" 




trimer (a) 







50 



100 



150 



path (AB) 



atom A (a) 



atom A (b) 




CPT 



trimer (a) 







50 100 

time 



150 



FIG. 1: (Color online) Heteronuclear trimer population for (a) 
S = —3 (in green) or (b) 8 — 3 (in red) and their corresponding 
numbers of atom A. Time is in units of A ; _1 , and 7=1. The 
dimer population (A2 or AB) remains essentially zero at all 
times and is not shown. The line labelled "CPT" shows the 
ideal, analytically derived trimer population. The evolution 
of the atoms B is also not shown. 

Figure 1 shows a numerical simulation of the formation 
of the heteronuclear trimers A2B, when considering the 
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channels AA and AB separately. In this specific example 
the atom A is 2 Na, with the s-wave scattering length 
3.4nm, and atom B is 87 Rb, with an s-wave scattering 
length of 5.77nm, A ; =4.718 x 10 4 s" 1 (I = 1, 2), and 

nj(f) = ni, sech(t/T), (4) 

with S\o/A; = 20, A;r=20. The s-wave scattering length 
for Na-Rb collisions, which depends on the details of the 
interatomic potential, is not well known. Here we take 
Xa=0.3125, Xb=0.5303, Xab=0.4214 and the other colli- 
sions parameters are taken as 0.0938 (in units of Xi/n) 
[20] . Very little is known about the scattering lengths of 
collisions involving molecular trimers, so we have carried 
out simulations with several sets of plausible parameters 
[20]. We found that the stable formation of trimers is 
always possible for a range of values of the detuning 5. 

The curves labelled (a) in Fig. 1 give one such example 
for <5 = — 3. Their general features resemble those of Fig. 
2 in Rcf. [4], which corresponds however to the creation 
of dimers rather than trimers. In particular, we observe 
a similarly increasing departure of the population trans- 
fer from the CPT solution as time evolves. Just as is 
the case for dimer formation [3] , the association of atoms 
into trimers is characterized by the existence of regions 
with unstable dynamics. This is for instance the case 
for for (5 = 3 (curves (b) in Fig. 1.) As expected, the 
two reaction channels lead to different dynamical behav- 
iors. In particular, the AB channel leaves a significantly 
larger number of atoms A in the sample at the end of the 
conversion process, and hence results in a lower yield of 
heteronuclear trimers. 

We conclude the discussion of the single-channel cases 
by noting that we have also studied the full quantum dy- 
namics of heteronuclear or homonuclear trimcr formation 
in the short-time limit using a c-number positive- P rep- 
resentation approach [6]. We find that in that limit the 
quantum dynamics reproduces the CPT dimer produc- 
tion predicted by the mean- field theory [21], and quan- 
tum noise-induced trimer damping occurs only near a 
total atom-trimer conversion [6, 21]. 

We now proceed to demonstrate that when acting in 
concert, the two channels can yield a significantly larger 
conversion rate and approach the ideal CPT yield of 1/3, 
see Fig. 2. Here, the coexistence of the two channels 
provides considerable additional flexibility in approach- 
ing the ideal CPT value for trimer formation. Note how- 
ever that this approach either requires an accidental co- 
incidence of FRs for the A2 and AB dimer formation, 
or might be realizable in other cases by applying a mag- 
netic field gradient [22] across the coexisting A and B 
condensates. 

The mean-field description of the two-channel situation 
is a straightforward extension of the single-channel case. 
Requiring as before that the number of dimers remains 
equal to zero, the CPT steady-state number of trimers is 
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FIG. 2: (Color online) Normalized CPT trimer population as 
a function of r\i = X1/Q1, 1 = 1,2. Also shown are three values 
of R — 772/171. At constant ratios fi2 (f)/Oi (i) the counter- 
intuitive evolution of the system is along lines of constant R 
starting from the origin. 

then 

N = (x 1 /n 1 )(\ 2 /n 2 ) 2 

9 " A 1 /0 1 + A 2 /fi 2 + 3(A 1 M)(A 2 M) 2 ' [) 

where the asymmetry between the two channels results 
from the fact that the intermediate dimer involves two 
undistinguishablc particles in the first one and two dis- 
tinguishable particles in the second case. 

Figure 2 plots the steady-state trimer number N g>s as 
a function of r\i = Xi/Cti, I = 1, 2, as well as the parameter 

R = m/m- (6) 

(Note that there is no CPT solution for R < 0.) As the 
STIRAP PA pulses f2 x (i) and Q 2 (i) are applied, 771 and 
772 increase and if the ratio of their amplitudes remains 
constant the system evolves along a line of constant R. 
To determine whether an optimum choice of R permits 
to approach the ideal trimer population of 1 /3 under the 
non-ideal STIRAP conditions of our nonlinear system we 
have solved numerically the mean-field equations of mo- 
tion of the system for various values of R £ [1, 3] (see Fig. 
3), using the same parameters as in the single-channel 
case. 

We found numerically that R = 2 leads to a trimer 
production that most closely approaches the ideal CPT 
solution, and is significantly larger than in the single- 
channel situation of Fig. 1. This is illustrated in Fig. 3c 
for (5 = and 5 = — 3. Note the insensitivity of trimer 
production to the detuning in that case. These results 
should be contrasted to Fig. 3b, which shows the evo- 
lution of the trimer population in case R = 1, again for 
(5 = and S = —3. Here, the trimer production is very 
significantly reduced, and depends strongly on the value 
of the detuning. Similar results have been obtained for 
the other values of R and the other sets of collisions pa- 
rameters [20] that we have considered. 
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FIG. 3: (Color online) (a) Two-channel generation of het- 
eronuclear trimers. The time dependence of the populations 
of the trimers (solid) and of atom A (dotted) are shown for 
6 = -3 and 8 = and (b) R = 1 or (c) 2 . The CPT value of 
trimers is also plotted. Xd x d 2 = and di,2 denotes A2 or AB. 



We observe also that in the two-channel case the trimer 
population can reach a transient value that is larger than 
its final value. This suggests that maintaining a constant 
ratio R during the evolution of the system may not be 
optimal. Future work will use genetic algorithms to de- 
termine the optimum time dependence of R(t) for maxi- 
mum trimer production. 

In conclusion, we have shown that a STIRAP scheme 
based on Feshbach-assistcd photoassociation, which has 
previously been shown to result in the production of ul- 
tracold molecular dimers [3] , can be extended to the gen- 
eration of molecular trimers, and that in the case of het- 
eronuclear molecules the interference between two for- 
mation channels can lead to a significant enhancement 
of trimer production. Future work will study further 
ways to optimize this scheme, and will give a quantum- 
mechanical description of the system, in particular with 
the goal of understanding the role of quantum fluctua- 
tions in the early stages of trimer production and the 



quantum statistics of the trimer field. We also plan to 
study the stability properties and the adiabatic geomet- 
ric phase of the system [4] , the superchemistry reaction 
A2+B — ► AB+B via a dark state technique, and possibly 
trimer creation or amplification in an optical lattice [23] . 
While experiments along the lines of this analysis promise 
to be challenging, recent progress in Efimov superchem- 
istry [9-12, 15] and in the manipulation of dimer-atom 
resonances [11-14] indicates that they may become pos- 
sible in the not too distant future. 

This research is supported in part by the US Office of 
Naval Research, by the US National Science Foundation, 
by the US Army Research Office, by the Joint Services 
Optics Program, by the National Aeronautics and Space 
Administration and by the NSC and SCUT. The authors 
thank Dr. M. Bhattacharya for helpful discussions, and 
Dr. P. Julienne for bringing Ref. [7] to their attention. 
* Electronic address: pierre.meystre@optics.arizona.edu 



[1 
[2 

[3; 

[1 

[5 
[6 

[r 
[8 

[9 
[10 
[11 

[12 

[13 
[14 
[15 

[16 

fi- 
lls 

[19 
[20 



[21] 



J. R. Anglin and W. Ketterle, Nature (London) 416, 211 
(2002); P. Meystre, Atom Optics (Springer- Verlag, 2001). 
T. Kohler, K. Coral, and P. S. Julienne, Rev. Mod. Phys. 
78, 1311 (2006); K. M. Jones, el al., ibid. 78, 483 (2006). 
H.-Y. Ling, H. Pu, and B. Seaman, Phys. Rev. Lett. 93, 
250403 (2004). 

H. Pu, et al., Phys. Rev. Lett. 98, 050406 (2007). 

K. Winkler, et al., Phys. Rev. Lett. 95, 063202 (2005). 

J. J. Hope and M. K. Olsen, Phys. Rev. Lett. 86, 3220 

(2001). 

R. Amado and J. Noble, Phys. Rev. D 5, 1992 (1972). 
V. Efimov, Phys. Lett. B 33, 563 (1970); Sov. J. Nucl. 
Phys. 12, 589 (1971); Com. Nucl. Phys. 19, 271 (1990). 
T. Kraemer, et al., Nature (London) 440, 315 (2006). 
M. Y. Yamashita, et al, Phys. Rev. A 68, 033406 (2003). 
P. Soldan, et al., Phys. Rev. Lett. 89, 153201 (2002); M. 
T. Cvitas, et al., ibid. 94, 200402 (2005). 
E. Nielsen, H. Suno and B. D. Esry, Phys. Rev. A 66, 
012705 (2002). 

M. Mackie, et al., Eur. Phys. J. D 31, 273 (2004). 

C. Ospelkaus, et al., Phys. Rev. Lett. 97, 120402 (2006). 

E. Braaten, H.-W. Hammer, and M. Kusunoki, Phys. 
Rev. Lett. 90, 170402 (2003). 

P. D. Lett, et al., Phys. Rev. Lett. 71, 2200 (1993). 

F. Colavecchia, et al., J. Chem. Phys. 118, 5484 (2003). 
K.-P. Marzlin, W. Zhang, and E. M. Wright, Phys. Rev. 
Lett. 79, 4728 (1997). 

J. Cheng, et al., Phys. Rev. A 73, 035601 (2006). 
For example, in the Thomas-Fermi approximation, we 
can take Na-Rb scattering length as half the critical value 
[(M Na M Rb ) 1/2 /(M N a + Af Rb )](a Na a Rb ) 1/2 for the stable 
coexistence of Na and Rb BECs. Here MN a ,Rb and ONa,Rb 
denote the mass and scattering length of Na and Rb 
atoms, respectively. Thus, by using Xab~0.0832 and as- 
suming the other collisions parameters in both cases to 
be those of atom A, except for x<T5=- 1-4583 (o~=a, b and 
q=d, g), we can reach a similar result as Fig. 1 (with 
fewer numbers of trimers). 

Details of the derivations will be published elsewhere. 



[22] H. Xiong, et al, Phys. Rev. Lett. 95, 120401 (2005). (2004). 
[23] F. Zhou and G. W. Semenoff, Phys. Rev. Lett. 97, 180411 
(2006); C. P Search and P. Meystre, ibid. 93, 140405 



